Green's Function Method for Line Defects and Gapless Modes in Topological 
Insulators : Beyond Semi-classical Approach 



O 
(N 

m 



Ken Shiozaki and Satoshi Fujimoto 

Department of Physics, Kyoto University, Kyoto 606-8502, Japan 
(Dated: January 16, 2012) 

Defects which appear in heterostructure junctions involving topological insulators are sources of 
gapless modes governing the low energy properties of the systems, as recently elucidated by Teo and 
Kane [Physical Review B82, 115120 (2010)]. A standard approach for the calculation of topological 
invariants associated with defects is to deal with the spatial inhomogeneity raised by defects within a 
semiclassical approximation. In this paper, we propose a full quantum formulation for the topological 
invariants characterizing line defects in three-dimensional insulators with no symmetry by using the 
Green's function method. On the basis of the full quantum treatment, we demonstrate the existence 
of a nontrivial topological invariant in the topological insulator-ferromagnet tri-junction systems, 
for which a semiclassical approximation fails to describe the topological phase. Also, our approach 
enables us to study effects of electron-electron interactions and impurity scattering on topological 
insulators with spatial inhomogeneity which gives rise to the Axion electrodynamics responses. 

PACS numbers: 73.20.-r, 73.43.-f 
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I. INTRODUCTION 

Topological phases realized in condensed matter sys- 
tems have been currently attracting much interest 
The topological classification of the ground states of 
bulk band insulators or superconductors is a useful ap- 
proach for the elucidation of low energy properties of 
the topological phases. The fundamental consequence 
of this classification is the bulk-boundary correspondence; 
i.e. there exist the topologically protected gapless 
states on a boundary surface which separates two sys- 
tems belonging to different topological classes i^"^'^^ Re- 
cently, Teo and Kane proposed the generalization of this 
topological classification to arbitrary topological defects 
such as line defects and point defects in insulators and 
superconductors 1^ It was suggested that the existence of 
gapless modes in defects follows from topologically non- 
trivial ground states of Hamiltonians which vary with 
material parameters characterizing the defects. The Teo- 
Kane theory was successfully applied to various kinds 
of heterostructure junction systems, describing correctly 
low energy gapless modes which appear at junctions. A 
basic idea of Teo-Kanc theory is an " adiabatic deforma- 
tion" . The interface structure with a finite energy gap 
between the ground state and the first excited state can 
be deformed smoothly into the Hamiltonian varies slowly 
in real space without closing the energy gap. Then, deal- 
ing with the spatial coordinates r which parametrize 
defects as adiabatic parameters, one can introduce the 
adiabatic Hamiltonian H{k,r), which can be used for 
the argument on topological properties of defects. Al- 
though the underlying idea is quite general, and indepen- 
dent of approximation schemes of calculations, a simple 
method for implementing this idea is to apply the semi- 
classical approximation to such a slowly varying Hamil- 
tonian, and to treat the spatial coordinate r and the 
momentum k as independent variables. In fact, this ap- 
proach was adopted for the argument on the topological 



insulator-antiferromagnet junction in refi^. The semi- 
classical approach is particularly useful when the spatial 
variation of material parameters raised by defects is suffi- 
ciently slow. In the case of line defects with broken time 
reversal symmetry, which may be raised by the topologi- 
cal insulator-magnet junctions, the topological invariant 
within a semi-classical approximation can be represented 
by the winding number of the axion 0-terroi^ which char- 
acterizes the topological structures of band insulators in 
bulk regimes. As long as the spatial variation of material 
parameters in the vicinity of defects is sufficiently slow, 
we can calculate the 0-term by using the semi-classical 
approach , 15,39,40 

However, in some cases, a semi-classical approxima- 
tion fails to give the correct energy spectrum of spatially 
inhomogeneous systems. For instance, in the case of a 
topological insulator-ferromagnet heterostructure junc- 
tion, which is particularly important for the application 
to the quantum Hall effect, a semi-classical approxima- 
tion which neglects quantum corrections leads to the clos- 
ing of an energy gap at the interface between the topolog- 
ical insulator and the ferromagnet. Since the full quan- 
tum energy spectrum of this system has an energy gap, 
the gap-closing is a spurious effect raised by the semi- 
classical approximation. Because of the gap-closing, we 
can not evaluate the topological invariant, and hence, the 
semiclassical approach for the topological classification of 
defects is flawed. 

Motivated by the above consideration, we propose a 
full quantum formulation for the topological invariant 
characterizing line defects in three-dimensional insula- 
tors with no symmetry exploiting the Green's function 
method, but, without using both adiabatic deformation 
and semi-classical approximation. On the basis of the 
full quantum treatment, we demonstrate the existence 
of a nontrivial topological invariant in the topological 
insulator-ferromagnet tri-junction systems of which the 
energy gap is closed within the semi-classical approxima- 
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tion. 

The organization of this paper is as follows: in SecHIl 
we briefly review the Teo-Kane classification scheme, par- 
ticularly focusing on the topological invariant for line de- 
fects in a three-dimensional insulator. We also demon- 
strate that in a topological insulator-ferromagnet het- 
erostructure system a semi-classical approximation ne- 
glecting quantum corrections to the energy spectrum fails 
to capture the correct topological features of the sys- 
tem. In Sec lIIIi we present the full quantum construc- 
tion of the topological invariant for line defects with 
broken time-reversal symmetry. In Sec lIVI and SecfVl 
we apply our approach to the topological insulator- 
ferromagnet heterostructure system. For this purpose, 
we use an exactly solvable model of the heterostructure. 
We present the exact solution of the heterostructure sys- 
tem in Sec lIVi and then, show the numerical results of 
the topological invariants in Sec|Vl In Sec lVIi we give a 
conclusion and discussions. 



II. THE TEO-KANE'S THEORY AND THE 
SEMI-CLASSICAL APPROACH 



the context of the axion electrodynamics describing the 
electromagnetic property of the topological insulatorsJ^ 
The above argument is based on a semi-classical ap- 
proach because quantum corrections to the energy spec- 
trum from spatial inhomogeneities are neglected. Hence 
it is not guaranteed that eq.(l) can describe the existence 
of gapless modes in cases where the semi-classical approx- 
imation fails. In fact, this may be the case when the scale 
of the spatial variation in the vicinity of interfaces can 
not be neglected. 
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The Teo-Kane theory assumes the semi-classical 
Hamiltonians H{k, r) which vary slowly in real space far 
away from defects. The existence of gapless modes local- 
ized at the defects follows from topologically nontrivial 
ground states of Hamiltonians H{k,s) := H {k,r{s)), 
where s parametrizes a line or a surface surrounding the 
defects. A similar semi-classical approach to the clas- 
sification of topological defects was also considered for 
the superfluid '^Hei^ In the cases of line defects in three 
dimensional insulators with no symmetry, the classifica- 
tion of the ground states of the semi-classical Hamiltoni- 
kz,s) is an integer Z characterized by the 
second Chern number 



Ch2 = --^ / tr[^A.F], 



87r2 



(1) 



here T = dA -I- ^ A ^ is the Berry curvature form^ as- 
sociated with the non-Abelian Berry connection Aij = 
{ui\duj) characterizing the valence band eigenstates 
|ui(fc,s)). The second Chern number C/12 can be rep- 
resented as the winding number of the O-term describing 
the local band structure 



2tt Jgi ds 



e{s)^2TT / d'^k Qz{k,s) 



(2) 



(3) 



where Q-i{k,s) is the Chern-Simons 3 form Q3 = 
- tr [A^dA+lA^A^A\. Then the nonzero 
winding number of 6 corresponds to the existence of gap- 
less modes. The 0-term is the same as the axion field in 



FIG. 1: (Color online) Tlie lieterostructure geometry for the 
topological insulator-ferromagnetic insuIator(FI) tri-junction. 

For instance, the semiclassical approximation fails for 
a topological insulator-ferromagnet heterostructure junc- 
tion. In the following, we demonstrate this by using the 
three dimensional Dirac model for the surface of the topo- 
logical insulator with ferromagnetic perturbation, 

H ^ v^ii (fca-CTi - idya2 - idlers) + m{z)fi3 + h{y)a3. 

(4) 

Here cr = (cri,CT2,(73) and /x = {fix, ^2, ^J■3) are the Pauli 
matrices representing spin, and orbital degrees of free- 
dom, respectively. The mass parameter 771(2) describes 
the difference between trivial and topological phases. 
h(jj) describe the domain wall of the ferromagnetic in- 
sulator. We consider surface geometry in which to > 
for z > , m < for z < 0, /i > for y > 0, and h < 
for y < as shown in FIG. [TJ The system has transla- 
tional symmetry in the a;-direction, and thus, is the 
only parameter of the Hamiltonian. This Hamiltonian 
possesses a chiral edge mode localized in the line defect 
at y, z ~ 0. The associated wave function is 



/ 1 



' dz m{z') /v — dz' h{y) / V ikxX 



(5) 



up to the normalization factor. Here we take 
the basis of spin and orbital degrees of freedom. 
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{1(73 = 1,M3 = 1), = l,Ai3 = -1), ks = -1,M3 = 1), 

I (73 = — l,/i3 = —1)}- Within the semi-classical approx- 
imation, p ~ — iV and x arc treated as commutative 
variables. The semi-classical Hamiltonian describing the 
interface of topological insulator-ferromagnet is 



H[kx, ky,kz', y, z) 

= Vfii {k^ai + ky(72 + fc^CTs) + m{z)^i3 + h{y)a3. 

The square of the energy eigenvalues are 

E (^kx, ky, kz, y , z) 



(6) 



v^{kl + kl) + ( ^m^{z)+v^kl ± h{y) 



(7) 



Near the interface between topological insulator and fer- 
romagnet in which the condition |m(z)| < is sat- 

isfied, there are the gap closing points at {kx,ky,kz) = 
(0, 0, ±-\/ft,2(y) — •m?{z)/v). Therefor the second Chern 
number ^ calculated from this incorrect spectrum can- 
not be quantized. This gap closing is a spurious one 
raised by the inaccurate semi-classical approximation. To 
obtain the 0-term for this junction system, we need quan- 
tum corrections which lead to the correct energy spec- 
trum. In the next section, we construct the 0-term in the 
full quantum manner. 

As discussed in ref.dHiil, in the case that a 
time-reversal symmetry breaking term is a Dirac type 
hafiy, z)^2 expressed by 7-matrices : 



H 



-iv^icr ■ V + m{y, z)/i3 -f- haf{y, z)^i2, (8) 



the semi-classical approximation is successful, since the 
energy gap does not close. The five set of ma- 
trices {/iiCTi, /xi(T2, A'lcrs, /^3, ^'2} consist Dirac matrices 
70: 71 1 72, 73, and 75. Due to the anti-commutative prop- 
erty, the nonzero m{y,z) and haf{y,z) play a role of 
a gap function, and set a characteristic length scale 



^ v/ ^^m'^(y, z) + h'^j:{y, z). To deal with the quan- 
tum corrections for the spatial inhomogeneity of m(y, z) 
and haf{y,z), the gradient expansion can be applied as 
long as the spatial inhomogeneities are slower than ^c- 
Especially, the semi-classical approximation, which is the 
zeroth order of the gradient expansion, is also applicable. 
In fact, the semi-classical energy spectrum of the Hamil- 
tonian (IHl), 



E[kx, ky,kz, y, z) 



±^v^{kl + kl + kl)+ m^{y, z) + hl^iy, z), 



(9) 



has an energy gap unless m?{y,z) -f h'^j:{y,z) = 0, i.e. 
away from the defect line, and hence, the second Chern 
number ([1]) is well defined. 

On the other hand, in the case of the topological 
insulator-ferromagnet junction, the Zeeman term ha^, 
which is not expressed in terms of 7 matrices, does 



not protect the energy gap since the ha^ commute with 
fczMi'''3- Thus the gradient expansion and semi-classical 
approximation cannot be applied to the interface of topo- 
logical insulator- ferromagnet. 



III. THE CONSTRUCTION OF THE 
TOPOLOGICAL INVARIANT FOR LINE 
DEFECTS WITH BROKEN TIME-REVERSAL 
SYMMETRY 

In this section, we present the full quantum treat- 
ment of the topological number characterizing line de- 
fects without time reversal symmetry in three dimen- 
sional Insulators. This can be achieved by using the 
Green's function. In spatially inhomo^eneous systems, 
the single particle Green's functions Gai3{i^^,Xi,X2) is 
defined by the Dyson's equation. 



dx2K{iuJ,Xi,X2)G{iuj,X2,X3) = (5(a;i — aJa). (10) 



K{iuj,Xi,X2) is the kernel of the Green's func- 
tion. For noninteracting systems K{iuj, Xi,X2) = 
[iuj — H{~i'Vi,Xi)] 5{xi — X2). In this cases, Dyson's 
equation is 

[iuj - H{-i'Vi,Xi)] G{iuj, Xi,X2) = 5{xi - X2). (11) 

Then we consider the Wigner transformation of 
G{iuj, Xi,X2)., 



G{iLu,p,R) := drG 



-tpr 



(12) 



Here, p is the momentum of relative coordinate r = 
Xi — X2, and R = {xi + X2)/2 is the center of mass 
coordinate. For translational symmetric systems, G is 
the matrix inverse oi K , G ~ K^^ , but in the case of in- 
homogeneous systems, G 7^ K~^. We introduce a closed 
path C — R{s) surrounding the line defect in real space. 
If the Green's function G{iuj,p, R) is not singular for all 
cj, p on the closed path C, then topological invariant N 
can be defined by 



N-- 



ds I du I d^p e'^'P"" 



tr [G-'^df,GG''^d^GG''^dpGG-^daGG~^d^G] , 

(13) 

where {n,v,p,a,r]) run over {ui,Px,Py,Pz, s), e is fully 
anti-symmetric tensor with e"P=:P«P^'' = 1^ and dg = 
^^^s^ • ^ is the directional derivative along the path 
C . This topological number characterize the homotopy 
of the map G : {uj,Px,Py,Pz, s) G{iuj,Px,Py,Pz,s) G 
GL(n,C), here n is the number of the band. If the base 
space {u!,Px,Py,Pz-, s) can be considered as this map 
is classified by the homotopy group tts (GL(n, C)) = Z, 
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which imply the existence of the nontrivial topological 
line defects. The same topological invariant as was 
previously obtained by Silaev and Volovik^ by using gra- 
dient expansions for the semi-classical Green's function. 
We stress that, in our approach, we have constructed 
the topological invariant (|13p from the argument on the 
homotopy of the exact Wigner-transformed Green's func- 
tion, without using gradient expansion, and thus Eq. (|13p 
is a full-quantum expression free from semi-classical ap- 
proximations. 

The topological invariant is not changed by the de- 
formation of the closed path C as long as the singularity 
does not occur on the path C. The nontrivial topological 
invariant iV 7^ means the existence of the singularity 
of G inside C, say, on the line defect. The singularity of 
Green's function G arises from the existence of ill-defined 
points in the map G, which correspond to or cxj eigen- 
values of the matrix G{iu!,p, R) at some points {i-u,p) on 
the line defect coordinate R. The 00 eigenvalues imply 
the existence of the gapless modes localized at the line 
defects because only w = conforms to the 00 eigenvalue 
for the definition of G. Also, as noted by 43ll44 [. in the 
case with electron correlation effects, the eigenvalue is 
allowed, and associated with the Mott insulating phase. 

The topological invariant N can be represented as the 
winding number of a certain potential function 6{R). 
This can be seen by rewriting ea. (jl3p as 



TV = — / dR{s) ■ A{R), 
27r Jc 



A{R) 



487r2 



tr [G-^dRGG-^d^GG-^d^GG~^dpGG^^d„G] . 

(14) 

According to the path-independent property of the line 
integral, A{R) is rotation free, V x A = 0. Hence A{R) 
can be represented as the gradient of a certain potential 
function e{R), A{R) = V6{R). Due to the Hermiticity 
of the Hamiltonian, A{R) is real, so 0{R) is also the real 
function. Thus we get 



N 



(15) 



In the cases of the slow limit of spatial variation 
or within the semi-classical approximation, the Green's 
function is reduced to a simple form, G{iuj,p,R) = 
[iu! — H{p, R)] . Then N is equal to the 2nd Chern 
number of the semi-classical Hamiltonian H{p, s) := 
Hip,R{s)): 



N = Ch2 



(16) 



in consistent with the Teo-Kane results.— 

In the semi-classical approximation, the above expres- 
sion of the 9 in terms of the Green's function is the same 
as the formula of the 6 term introduced by Qi, Hughes, 



and Zhang for the Axion electrodynamics of the topolog- 
ical insulators It is noted that eq. represents the 
topological invariant even when the semiclassical approx- 
imation, and the gradient expansion for the derivation of 
the Axion electrodynamics action are not justified, as in 
the case of the topological insulator-ferromagnet junction 
mentioned in Sec. II. However, since we do not know how 
to derive the Axion electrodynamics action without using 
the gradient expansiouj^i^ we could not make a direct 
connection between the 9 introduced above and the 9- 
term of the Axion electrodynamics in the case that the 
semiclassical approximation is not applicable. 

Here, we make a general remark on the calculations of 
topological invariant associated with defects. To avoid 
the singularity at the defects, the closed line or surface 
surrounding the defects must be enough far away from 
the defects. For example, in the case of line defects in 
heterostructure systems with broken time reversal sym- 
metry, we calculate the changes of the ^-term on the in- 
terfaces of insulators which are located far away from 
line defects by using eq.(ll). This procedure is the same 
as the calculation of the bulk topological invariants for 
which effects of boundaries are neglected. 



IV. AN EXACTLY SOLVABLE MODEL FOR A 
TOPOLOGICAL INSULATOR-FERROMAGNET 
JUNCTION 

In this section and the next section, we demonstrate 
that the formula ([T^ gives the correct topological in- 
variant even for the case where the semi-classical ap- 
proach fails. For this purpose, we use an exactly solvable 
model of a topological insulator-ferromagnet heterostruc- 
ture junction, which allows us the full quantum treat- 
ment of the spatially inhomogeneous system. We will 
calculate the variation of 9 on the surface of this system. 
The model Hamiltonian has a kink structure in the mass 
parameter m, which distinguish between topological and 
trivial phases, 

H = vfii {Kai + kyij2 - idzas) + m{z)fi3 + has. (17) 

We consider surface geometry in which m > for z > 0, 
and m < for z < 0. First we investigate the energy 
spectrum of this Hamiltonian. According to the appendix 
1X1 the energy spectrum of this Hamiltonian (flT)) can be 
labeled by the using of the eigenvalues e of the Jackiw- 
Rebbi Hamiltonian, H^iz) := —i^idz + m{z)fi3. The 
results are 



E,{k,, ky) = ± ^{e + h)^+v^{kl + kl), 
±^{e-hY+v^{kl + kl), 



(18) 



for e > 0, and 



Eo{k,,ky) - ±^v^k^^ + k^) + h^, (19) 
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for e = 0. As shown in Appendix in the case of 
m{z) = TOtanh(z/^), we can obtain the exact Green 
function as well as the exact energy spectrum.— The 
possible values of e are classified into the bound states 
localized at z ~ 0, and the continuum states e'^ 



0, 



(20) 



where v = ^ describing the spatial scale of the kink. 
When V > 1, bound states with nonzero energies are 
allowed as well as a zero energy state, which corresponds 
to = 0. The energy dispersion of the bound states 



E, 



{k,,ky) = ±y^(e„±-/i)2+v2(fc.2+fc2), 



(n = 1,2,-- - < I/) 



(21) 

implies that the energy gap is closed only when |e„-|-| = 
\h\. This property is consistent with the gap closing 
found in the semi-classical approximation because of the 
following reason. The semi-classical approximation cor- 
responds to the infinitely slow varying kink potential : 
^ = oo. Then the bound states distribute in the range of 
— m < e„ < TO densely, and hence the energy gap is closed 
because there are the bound states with |e„±| ~ \h\. This 
behavior is drastically changed in the cases of the topo- 
logical insulator-antiferromagnet heterostructure consid- 
ered in ref . [s^ . According to the Appendix El in the 
case of the anti-ferromagnet, the energy dispersion for 
the bound states are 



<±(fcx,fc,) 



± 



^l± + hlf + v-'{kl + kl), 
(n = l,2,.. 



(22) 



The system is the fully gapped with respect to the arbi- 
trary large ^, so the semi-classical approximation is valid. 

In real materials, ^ is atomic scale, and continuum ap- 
proximation is not valid for describing the interface struc- 
tures. However, we are mainly concerned with how the 
9 term changes at the interfaces. This change is a topo- 
logical property, and does not depend on the details of 
the microscopic structure of the interfaces as long as the 
energy gap is not closed. Thus, for simplicity, in the fol- 
lowing, we set v — 1 ^ = ii/to, for which the expression 
of the exact Green's function is considerably simplified. 

Here, we summarize the exact results for the Green 
function of the model ([T7)) . The detail of the deriva- 
tion is presented in Appendixes [Bl [Cl and |D] It is con- 
venient to introduce the polar coordinates with respect 
to {uj,vkx,vky), Lo — psm6,vkx = p cos 9 cos (j),vky ~ 
pcos9 sincf). Then, the Green's function is decomposed 



as 



G{p,9,(f);zi,Z2) 

U\^)V\9)W^Gip, zi,z2)W(J2V{9)(T2U{dp), 



(23) 



where 



[/(</>) 



V{9) 



V2 



(l+im), (24) 



and G{p,zi,Z2) is the Green's function for the one- 
dimensional Hamiltonian H{z) : 



H{z) = —ivdzP-iCTs + m{z)ii2 + ha^. 



(25) 



H{z) is diagonal in the spin space, H{z) ~ H^{z) © 
H^{z). Then the problem to calculate the Green's func- 
tion reduces to solving one-dimensional Hamiltonians in 
the orbital space. 



{z) ~ —ivdzP-i + TOtanh 
H^{z) = ivdzP-i + TOtanh 



Ai2 



Ai2 + h, 
h. 



(26) 



lih — Q these systems are equivalent to the Jackiw-Rebbi 
model^ in the one-dimensional chiral symmetric system. 
The Jackiw-Rebbi model possesses a zero energy soliton 
solution, and the finite h induces the shift of the energy. 
The soliton modes in the two sectors H'^iz) and H^{z) 
correspond to the helical Dirac fermions on the surface 
of the topological insulators. 

When V = 1 the eigenfunctions for 0-3 = t sector take 
simple analytical forms: 




sech 



1 



(27) 



i'x/ v^P^ -\-m^ J 



with eigenvalues eg ~ h, and ep,± — ±-\/ri 
respectively. Due to the completeness property of the 

set of eigenfunctions I'I'o' *^p.±}' Green's function 

&{p, zi, Z2) can be expressed by the spectral represen- 
tation as follows. 



G^{ip, zi,Z2) 

= E 



dp 
2^ 



<±(^i) 



*I±(^2) 



ip - eo 



(28) 



The Green's function of (T3 = ^ sector is obtained in a 
similar way. The Wigner transformation of G^, and G^ 
is given by, 
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G\ip,p,R)^ 
G^{ip,p,R)^ 



ip — h 



vp — im{p, —vp, R/v; h) 



{ip - hy - w2p2 _ ^2 \yp + im[p^ vp, R/v; h) {ip - h)/3{p, vp, R/v; h) 



1 



[ip + h)"^ — v'^p^ — im? 



{ip + h)(3{p, vp, R/v; —h) —p — irh{p, vp, R/v; —h) 
-vp + im{p,~vp, R/v; —h) ip + h 



(29) 



The derivation of and the definition of the effective 
mass TO, and the factor /3 are presented in Appendix [Dl 
After integrating over 9, and (j), we obtain (sec Ap- 
pendix |E] for details), 

MR.) = ^J^"" dp J dp. 



(30) 



where 



Xt,ip,Pz,Rz) ^ G ^df,G {pe {p,p.,R.)) 
Xe{p,Pz,Rz) = P-i(Ti +G^^piaiG , 
^^{p,p.,R.)=d^GG~' {pe{p,p.,R.)) 
^g{p,Pz,Rz) = piai+GpiaiG~^ , 



(31) 



G = G^^ S) G^, and e is fully anti-symmetric tensor with 
gPSp.i?. = 1 The derivation of Eq . ^ is in Appendix 
IeI Eqs. pop and ([3T|) are used for the numerical calcu- 
lation of 9 and the topological invariant in the next 
section. 



V. NUMERICAL RESULTS OF 6 AND THE 
TOPOLOGICAL INVARIANT FOR THE 
TOPOLOGICAL INSULATOR-FERROMAGNET 
JUNCTION 

In this section, we present the results of numerical cal- 
culations of the variation of 9, A9, in the vicinity of inter- 
faces of the topological insulator-ferromagnet tri-junction 
system, for which the semi-classical approximation is not 
applicable. The heterostructure system considered here 
consists of the topological insulator and the two ferro- 
magnets with the opposite direction of the magnetiza- 
tion. There are three interfaces; i.e. one between the 
topological insulator and the ferromagnet with the inter- 
nal magnetic field h > (interface I), another between 
the topological insulator and the ferromagnet with h < 
(interface II), and the other between the two ferromag- 
nets with h > and h < (interface III). The intersec- 
tion of these three interfaces forms a line defect. The 
Hamiltonian in the vicinity of the interface I or II is, 

II{kx, ky, z; /i) 

= vpi {kxcri + kya2 — idzcrz) + mtanh(z/^)^3 -t- ha-i- 

(32) 
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FIG. 2: (Color online) The numerical calculations of ylz(i?z) 
given by the Hamiltonian (|32p . Points in figure are the results 
of numerical calculation for sample points i?^. The smooth 
lines are given by interpolation of the points. Black (solid), 
red (dashed), and blue (dotted) lines correspond to ft = 0.1m, 
0.5m, and 0.9m, respectively. 



Here we set ^ = 1. First, we consider the interface 
I (/i > 0). In FigUl we show the numerical results of 
Az{R.) for this model calculated by using pQ]). (pO| . and 
([5T]) for three values of /i > 0: h = O.Ito, 0.5to, and 0.9m. 
Also, we show in Fig l3] the variations of 9 obtained from 
the integration of A^i?^), 1^9{R.) = J^^ A{R'.)dR'.. It 
is found that for each value of h, the total change of 9 
is TT. This results are consistent with the flbuik defined 
by the Chern-Simons 3 form. Here we introduce the no- 
tation "6'buik" to distinguish from the 9 defined as the 
potential of A. The ^buik describes the bulk structures 
of insulators. Since inversion symmetry is not broken in 
the bulk ferromagnet and the bulk topological insulator, 
the value of ^buik is fixed to O for trivial insulators, or tt 
for topological insulators. Our method provides a reli- 
able method for calculating the spatial change of 9 which 
connects the trivial value in the ferromagnet region and 
the nontrivial value tt in the topological insulator region. 

Next we consider the interface II {h < 0). Just at 
h = 0, the induced mass gap in the helical Dirac fermion 
on the surface of topological insulator disappears. Then 
the total change of 9 is allowed to be a different value 
from the case of /i > since the Green's function have 
singularities when the energy spectrum is closed. In fact, 
the relation, 



P3aiH{kx,ky, z; h)p^ai = H{~kx, ky, z; -h), (33) 
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FIG. 3: (Color online) The variations of 9 on the interface 
of topological insulator- ferromagnet heterostructure. Black 
(solid), red (dashed), and blue (dotted) lines correspond to 
h — 0.1m, 0.5m, and 0.9m, respectively. 



follows 



and 



G{iuj, 



-h), 



A,{R,-h) = -A,{R,--h). 



(34) 



(35) 



Hence the total change of for /i < is — tt. The sign of 
the change of 9 is determined by the sign of ft,, i.e. the 
sign of the induced mass gap 

Finally, we consider the variation of 6 on the domain 
wall of ferromagnet (the interface III), the Hamiltonian 
of which is given by 



= Vfj.1 {kx(Ji 



idya2 + k^as) + mfis + h{y)a3. 



(36) 



Here h{y) describes the domain wall structure of ferro- 
magnet. it is found that 6 docs not vary on the domain 
wall of ferromagnet since the relation, 



follows 



and 



^i3a3H{k^,,y,k^)^i3as = H{k,.,y, ~k^), 
M30-3G(jw, k^,py, fc^, Ry)pL3a3 

= G{iui, ka:,Py, -k:,,Ry), 



0. 



(37) 



(38) 



(39) 



The total change of 6 about the closed path C 
surrounding the line defect in topological insulator- 
ferromagnet tri-j unction is tt + + tt = 27r, then we get 
the topological invariant N = 1. The nontrivial topolog- 
ical invariant means the existence of the singularity of G 
on the line defect. Then our result is consistent with the 
existence of chiral edge mode described by eq. ([5]). 



VI. DISCUSSIONS AND CONCLUSION 

In this paper we proposed a full quantum formulation 
for the topological invariant characterizing line defects in 
three-dimensional insulators with no symmetry by using 
the Green's function method. The nontrivial topological 
invariant leads to the existence of gapless modes local- 
ized in the line defects. Our approach is applicable to 
various heterostructure systems involving topological in- 
sulators, including the case that a semiclassical approx- 
imation for spatial inhomogeneity raised by defects fails 
to describe topological features of the systems. As an 
example of the failure of a semiclassical approximation, 
we considered the topological insulator- ferromagnet junc- 
tion system, which is important for the application to 
the quantum Hall effect. To deal with the issue beyond 
semi-classical approximations, we employed an exactly 
solvable model of the heterostructure junction, for which 
the exact Green's function can be obtained. On the ba- 
sis of the Green's function representation of the topolog- 
ical invariant, we demonstrated the nontrivial topolog- 
ical invariant for the topological insulator-ferromagnet 
heterostructure junction, which implies the existence of 
a chiral edge mode leading to the quantum Hall effect. 
Our method provides a reliable method for calculating 
the spatial change of 6 which connects the trivial value 
in the ferromagnet region and the nontrivial value tt in 
the topological insulator region. 

Within the semi-classical approximation, our formula 
of the topological invariant for line defects and the asso- 
ciated 9 term expressed in terms of the Green's function 
are, respectively, the same as the Teo-Kane's formula and 
the Qi-Hughes-Zhang's formula for the 9 term of the Ax- 
ion electrodynamics. Since we do not know how to derive 
the Axion electrodynamics action without using the gra- 
dient expansion, we could not make a direct connection 
between the 9 introduced in this paper and the 0-tcrm of 
the Axion electrodynamics in the case that the semiclassi- 
cal approximation is not applicable. We believe that even 
in such a case, the 9 expressed in terms of the Green's 
function is related to electromagnetic response functions 
characterizing the Axion electrodynamics of the topolog- 
ical insulators. The clarification of this point should be 
addressed in the near future. 

There are several future directions of the current study. 
For example, our approach based on the Green's func- 
tion may be useful for the study on effects of electron- 
electron interactions and impurity scattering in the topo- 
logical insulator with spatial inhomogeneity. Another 
direction is the application to the case of carrier-doped 
topological insulators, in which bulk metallic states co- 
exist with surface gapless states, as realized in various 
three-dimensional topological insulators J^ii^ii^ Recently, 
the Axion response in such gapless systems has been 
discussedj^i^ For such a case, effects of impurity scat- 
tering on the Axion response involving the bulk gapless 
states may become important, and can be evaluated by 
using our Green's function method. 
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Appendix A: The Energy Spectrum of the Hamiltonian (|17p 

Hamiltonian ()17p is decomposed into two parts: 

H{h^,ky-z) = Hi{k^,ky) + H2{z) 

Hi{kx,ky) = v^i{kxC7i + kya2) + haz (Al) 
H2{z) = -ivnidza^ + m{z)fi3. 

H2{z) is diagonal in the spin space, H2{z) = hI{z)(BH2{z). Due to the symmetric properties of H2{z) : ^2-^2(^)^2 = 
—H2{z), and (■^) ~ M3^2 ('^)M3i the eigenfunctions of H2{z) with nonzero eigenvalue can be represented as a set of 
four eigenfunctions with energy e, — e, e, — e : 

u^lJ , (^>0), (A2) 



y ' \^ ) ' yiJ.'iXeiz) ) ' \^llx<^{z) 

here Xe{z) is the eigenfunction of Hliz) with the eigenvalue e : H2{z)xt{z) = ext{z). Then H{kx,ky; z) is closed in 
the subspace spanned by ^e{z), 

H{kx, ky, z)$£(z) = ^t{z) [vfii{kx<yi + kya2) + has + • (A3) 

By diagonalizing of v^i(kx(Ji -t- ky(T2) + ha^ + e/xs, we get the energy spectrum labeled by e, 

E,{k,,ky)^±^{e + hY + v^{kl + kl) , ±^{e-hY+v^{kl + kl). (A4) 
Furthermore, -^2(2) have the pair of zero energy bound states. 



/I 




(A5) 



Mz) 

' Voy 

H{kx, ky, z) is closed in the subspace spanned by $0(2), 

H[kx, ky,z)<^o{z) = $0(2:) [v[kxT2 - kyTi) + hTs] . (A6) 

Here the Pauli matrices t are defined in the subspace spanned by ^q(z). By diagonalizing of v^k^Ti + kyT2) + hr^, 
we get the energy spectrum for zero energy bound states of H2{z), 



Eo{kx,ky) = ±^v^{kl + kl) + h\ (A7) 

The eigenfunctions $o(^) corresponds to the degrees of freedom for helical Dirac fermion on the surface of topological 
insulator. 

Here we remark the cases of topological insulator-antiferromagnet heterostructure described by the Hamiltonian 
considered in ref . [soj . 

H^^{kx, ky] z) = v^iiik^ai + feyO-2) + haf^J.2 - ^v^Jil^z(7z + m{z)^3. (A8) 
Then the Hamiltonian H^^ [k^, ky \ z) is closed in the same subspace spanned by $e(z), 

H^^{kx,ky, z)$e(z) = <^e{z) [vfiiik^ai + kya2) + hafH2 + em] ■ (A9) 
The energy spectrum energy spectrum labeled by e is 

Ef{k.,, ky) = ±^e^ + hl^+v^kl + k^y). (AlO) 

Here the upper and lower bands are doubly degenerated, respectively. The energy spectrum for zero energy bound 
states is the same form for topological insulator-ferromagnet heterostructure. 
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(B2) 



Appendix B: The Polar Decomposition of Green's Function 

The kernel of the Green's function associated with Hamihonian ([TT]). 

K{ll:, kx, ky] z) ^iu: - vfii {k^ai + kya2 - idzC^a) - m{z)fi3 - ha^, (Bl) 

is decomposed in the following way. We introduce the polar coordinates with respect to (w, vk^, vky), uj = pcos 9, vk^ = 
p sva 9 cos 4>,vky ~ psm9 sincj). Then (|B1[) is decomposed as 

K[p, 9, (p; z) = ipcos9 — vp,i {psin9 cos4>ai + p sin sin 0(72 — idzcr^) — m(z)ii3 — ha^ 
= ?7^((/)) [ipcos6' — psm9piai + ivdzPio^ — m{z)p:} — hcr^] U{4>) 
= U\(f))a2 [ipcos9(72 + ipsm9pi(73 — vdzPi<7i — ■m{z)p.^ij2 — ihai] U{4>) 
= U^(l))a2V\9) [ipG2 - vdzl^iai - m{z)pz<j2 - iha^] V{9)U{<j)) 
= U^i^)a2V^i9)a2 [ip + ivd^Pia^ - 777,(2)^^3 - ha:^] V{9)U{4>) 
= U\(t>)a2V\9)a2W^ [ip + ivdzPia^i - m{z)p2 - haa] WV {9)U {(j)) , 

where 

U{(j)) = ei"'''', V{9) = e^^^''^\W = ^(1 + ip^). (B3) 

Green's function is also decomposed as 

G(p,0,0;zi,z2) - U\<i,)V\9)W''G[p,z^,z2)WG2V[9)G2U{<l>). (B4) 

G(p, 2:1,22) is the Green's function for the one-dimensional Hamiltonian H(z), 

H{z) = —ivdzPi<7z + m{z)p2 + ha^- (B5) 

H{z) is diagonal in the spin space, H{z) = H^{z) H^{z). Each of sectors forms the Jackiw-Rebbi model with the 
constant energy shift ±h. 

Appendix C: An Analytically Solvable Model for the Jackiw-Rebbi problem 

The Jackiw-Rebbi model with kink mass m{z) — 777 tanh(2/^) 

H = —ivdzPi + 777 tanh (2/^) fi2 

-^(wa^ -1-777 tanh (2/0)^ (CI) 

—7 (vdz — Tretanh (2/^)) 



is analytically solvable. In this appendix we briefly sketch how to solve Hamiltonian (jCl|) according to ref. |45ll46j . 
The eigenvalue equations for $(2) = * {u+{z), u-{z)) are 

f-7(7;a^ + 777 tanh(2/^))71_(2) = £77+ (2), ^^^-j 

I —7 {vdz — Tretanh (2/^)) 77_|_(2) ~ eu-{z). 

For e 7^ we obtain 

v^d^z +e^-m^ tanh2(2/0 T ^ (l - tanh^{z/0) u±{z) = 0. (C3) 



u±{X) - 0, (C4) 



By introducing a new variable X = tanh(2/^) , eq. ()C3P may be rewritten as 



X2 



10 



where v = — ,/^^ = ^{m^ — e^). The solutions of these differential equations are represented by Legendre functions 



u+{X) ^ P^,XX), 
u^{X)r^Pt^{X), 



(C5) 



Pi^{X) 



1 fi + x 



r(i VI 

where F is the hypergeometric function, 



F 



1-X 



(C6) 



F[a,P,j,z] = J2 



71=0 



(7)n n\ ' 



(C7) 



(a)„ = a(a + l)(a + 2) • • • (a + n - 1). 



The relative phase factor between u^(X) and U-{X) is determined by ea. (|C2p . The eigenstates of which fj, is pure 
imaginary correspond to continuum states, and can be labeled by the wave number p, and a sign of the energy ± : 



1 



(C8) 



/2 \ ±i- 



-F 



Here, the normalization of $p,±(z) is determined by the asymptotic behavior at z — > oo, 

V 2 V ^m2_|_^,2p2 J 

The bound states are determined by the expression of P}^{X) as follows, 



v{l-^i) 



1 - ^ + -/i - i^, 1 — ^; ■ 



This follows from the formula. 



F [a, /3, 7; ^] = (1 - zy-'^-Pp [7 - a, 7 - /3, 7; ^] ■ 



(C9) 



(CIO) 



(Cll) 



To ensure that F is finite for X — > — 1 (i.e. for z — > —00), we impose fi < 0. F must be a polynominal of a certain 
degree n, say, — /i — v = —n. Then the non zero energy bound states can be labeled by integer n — 1,2, ■ ■ ■ , < v : 



n ( n 
e„.± = ±m^^ - 2 



2 J ' 



-r^ (1 - X^) — F \l-n,2v~^n,l^v-n; ^\ 

T» /I _ y2^,;^ Y^ri-l (l-Tt)fc(2i/-w)fc 1 / 
y^2i./n-l'' ^^==0 (l+i^-n). fc! I 2 ^ 

n _ \"2\ii^ Y^n (-")fc(l+2t/-«)fc 1 i \-X \^ 

up to normalization. The zero energy bound state 3>o(-2) «: e"-'^" '"(^ ~ 1/ [cosh(z/^)]'' is 

£0 = 0, 

$o(^) 



(C12) 



o\ /r(!. + i/2) 



(C13) 



V V e7ri/2r(zy) [cosh(z/eF' 
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In the case of = 1 the expressions of eigenstates are considerably simphfied. The continuum states are 



= — ^ I _|_ vp+tm tanh( ^ j e* 



and the zero energy state is the only allowed bound state 

£0 = 0, 



0\ 1 (C15) 



1/ V 2^; coshfsi^) 



Appendix D: The Wigner Transformation of Green's Function G(p, ^1,22) 

In this appendix we perform the Winger transformation of G(p, zi,Z2). For simplicity we set v = 1. We calculate 
for the 0-3 = 1 sector. The calculation for the 0-3 = 1 sector can be done in a similar way. The Green's function (pS)) 
is 

G^{ip, zi, Z2) 

dp' ( ip~h p' - ZTOtanh(TOZ2) \ ei-p'(^i-^2) 

^ \p' + »77ltanh(mzi) (»p - fe) (p'+»'»tanh(>np))(p|-.mtanh(^.,)) I _ _ _ 







m 



sech(TOZi )sech(TOZ2 ) . 



The Wigner transformation of is defined by 

G\ip,p, R) := 1 dr G\ip, R+^,R- \)e-^^' . (D2) 

In the integral of r, the nontrivial contributions stem from r-depending terms in G^{ip, R+ ^, R — There are four 
terms as shown below: 

(A) [ dr [ — ^"^tanh [m{R + r/2)] ^^'r^-jpr 
J J 2n [ip — hY — p'"^ ~ m? ' 

iB) [ dr [ — ^^"^^ ^"^'■^ ~ '^^^^^ c"''''c~"'^ 

J 7 27r (ip — hY — p'^ — ra^ ' 

(C) [drf^ (tp-h)e'P'^e-'P' (D3) 



2'K [{ip — hY — p'^ — n^^] (p'^ -l~ rn^) ' 
• [imp' (tanh [m{R + r/2)] - tanh [m{R - r/2)]) + rrY tanh [m{R + r/2)] tanh [m{R - r/2)]] 

/777 
dr -sech \m{R + r/2)] sech \m{R - r/2)] e-'^"". 
2{ip - h) 

The p'-integral can be performed by using the following formulae, 4 

dp' e'P''' _ e-^f''! 
2^ p'2 + ^2 ~ 2A ' 

2^ (p'2 + A^){p'-^ + m2) " "2^^°^'''' {ip-h)^ ' ^ ^ 

dp' e'-P" 1 /e-™l'~l e-'^l''l 



27r (p'2 + ^2)(p'2 ^2) 2{ip-h)^\ m A 
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where A = m? + {p + ih)^. Then, the four terms in eq. (|D3[) are rewritten into the foUowing forms, 



(A) 
(C) 

(D) 



ITfl f 

-— / dr tanh [TO(i? + r/2)] e~ 



A\r\ —ipr 



2A 



dr tanh[m(i?-r/2)]< 



-yl|r| —ipr 



2{ip~h) 



dr sgn(r) {tanh [miR + r/2)] - tanh [m{R - r/2)]} (^e^"!"! - e^'^l'^l) 



f tanh [m(i? + r/2)] + tanh [m(i?- r/2)] ^ 
2{ip-h) J ^ \ tanh(2TOi?) 



m 



e 

-A\r\ 



ipr 



-ipr 



dr 



-ipr 



(ip — h) J cosh(2m/?) + cosh(r 



(D5) 



We use the relations tanh(a;) tanh(j/) = '^'^"t'inh(a*+^^^^'* ^^'^ cosh(a;) cosh(y) — [cosh(x + y) + cosh(x — y)] /2. After 
performing the Fourier transformation with respect to r, we obtain, 



(A) 
(B) 
(C) 



^2A 
i 

2A 



f mR, - / -mR, 



/ ( mR, ^_^) _ / (mR, 



1 



2{ip- h) 



fimR,^^]+f(-mR,'^^]+f(mR,^I^]+f(-mR ' 



ft o A + ip \ ( „ A-ip \ ( „ A~ip \ ( „ A-ip 
-f mR, - / -mR, - / mR, - / -mR, 



2(ip- /i)tanh(277ii?) 



m + ip\ ( m-ip\ I m.-ip\ I m - ip 

f mR, + / -mR, - / mR, - / -mR, 



^ (/ Lr, ^) + / (-mR, ^) - / Li?, ^) - f (-mM, ^ 
A I \ m J \ m J \ m J \ m 

m^(ip — h) 



[{ip — hY — p^ — m^] {p^ + m^) ' 

27rsin(2pi?) 
{ip — h) sinh(2mi?) sin {■np/m) 



(D6) 



Here /(a, h) is defined by 



f{a, h) := / ds tanh [a + "-] e'^' = (-1 + 2i^[l, 6, 6+1, -e-^"]) , (3?(6) > 0), 



(D7) 



where F[a, /?, 7, z] is the hypergeometric function. F[\,b, 6 + 1, z] is the analytic continuation of the infinite series, 



F[l,h,h+\,z] = 5^ 



6 + n 



(|z|<l). 



(D8) 
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We set an effective mass 'm{p,p, R; h) and an effective factor /3{p,p, R; h) as follows, 



rh{p,p,R; h) 
p{p,p,R-h). 



P^ + A^ 

2A 
p^ + 

2{ip-hf 
A + ip 



f mR, 



A + ip 



f mR. 



m 
m + ip 



f -mR, 



A — ip 



m 



f -mR, 



m 

m — ip 



m 



f mR, 



m — ip 



/ -mR, 



m — ip 



m 



f mR 



m 
p^ + 



2(ip- /i)2tanh(2mi?) 



/ y-mR, 
f (mR. 



m 
1 



f mR, 



A + ip 



+ f ~mR, 



A — ip 

m 
m + ip 

m 

A — ip 



f mR, 



A — ip 



m 



f [-mR, 
- f (mR. 



m 
m ~ ip 

m 
A — ip 
m 



f -mR, 



A — ip 



f mR, 



f -mR. 



m 
m ~ ip 
m 

A — ip 



f -mR, 



m — ip 
m 



m 



p — m" [ip — hy — p — m 



27r sin(2pi?) 



p'^ + 



{ip — hy sinli(2mi?) sin (irp/m) 



where A = \/m'^ + {p + ih)'^ which depends on h. The Green's function for as = t sector is represented as 



G^iip,P,R) 



ip — h p — ifh{p,—p,R;h) 

{ip - hy - p^ - m^ \p + im{p,p,R;h) {ip - h)l3{p,p, R;h) 



(D9) 



(DIO) 



The effective mass m and effective factor /3 describe the spatial inhomogencity near the kink structure. In a similar 
way the Green's function for (73 = J, sector is represented as 



G\tp,p,R) 



1 



f {ip + h)(3{p,p, R; -h) ~p - im{p,p, R; -h) 

{ip + hy - p2 _ ,7^2 \-p + im{p, -p, R; -h) ip + h 



(Dll) 



Appendix E: The Derivation of Ea.([30l) 

According to the polar decomposition the integral with respect to </>, and can be performed. First, we set 
G{p,e,4>,p„R,) = U\4>)G'{p,e,p„R,)U{4>), G'{p,9,p„R,) = V^{e)W^G{p,p„R,)Wa2V{e)a2. Then 



G~^d^,G = U\^)G'-^d^G'{p,e,p„R,)U{^) , {fi^p,e,p„R,) 
G-'d^G = ^[/t(0) [cT3 - G'-VgG'] U{q^) . 



(El) 



The integrand do not depend on 0, 

A,{R^) = J dpd(j)d9dp, e^'^'P" tr [G-^d^GG-^d^GG-^dyGG-^dpGG-^d^G] 

= / dpdddp, t^^'P" tr [{0-3 - G"-V3G"} G-^d^G'G'-^d^G'G'-^dpG'G'-^d^G'] 

" 4^ / '^P'^^'^P' ^^""^ t'^s {G-^d^,G'G'-^d,,G'G'-^dpG'G'-^d„G' 

-dpG'G'-^d^G'G'-^dpG'G'-^daG'G-^)] , 

where {p,, v, p, u) run over {p, 9,pz, Rz), and e is fully anti-symmetric tensor with ^p^P'^' = 1. Next, we calculate the the 
integral with respect to 6. Due to the relations, a2V{9)a2 = cr^ (cos | + ipiCTi sin |) (J2 = cos | — ipiCTi sin | = 1^^^), 
we obtain, 



G'-^dpG' ^ v{e)w^G-^dpGWV\e), {p e {p,Pz,Rz)) , 

(E3) 



G'-^deG' = --V{0)W^ l^piai + G^Vif^iGj WV\9) , 



dpG'G'-^ = v\e)wdpGG-^wW{e), {p e {p,pz,Rz)) 



deG'G'^^ = --V\e)W \piai + Gpi<7iG-^\^ W/V(6') 
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Here, we introduce Xtn^i-i to simplify notations as follows, 



Xeip,Pz,Rz) = ^J■l(Tl+G ^HiaiG , 
^t,{p,Pz,Rz) ^ df,GG , {iJ.&{p-,Pz,Rz)) , 



Then, we have. 



Az{Rz) = j dpdOdp, ^^-P" tr h {V(B)x^.XuXpXay\Q) - V\Q)i>.i.ipiaV ■ (E5) 

For V^6)a3V{6) = (cos | — i^iCTi sin |) 0-3 (cos | + i^iCTi sin |) = 0-3 (cos6' + i/.Ji(Ti sinfl), and V{9)a3V^9) = 
(cos I + i/ii(Ti sin |) (73 (cos | — ipi<Ji sin |) = (T3 (cos0 — ipi<7i sinfl), eq. (IE5[) is rewritten into, 

^z(i?2) = -^ J dpdp^ J de e'^'^P'^ tT[a3{{cos9 + ipiaisme)xtj.XuXpX<j-- {cos9-ifiiaisme)(^^,,^p^„}] 

= y rfprfp. e^"^" tr [(73 {2i//i(TiXMXi^XpX- + Sz/iKTi^^e.epe^}] (E6) 

= ^ J dpdp^ e'^'^P'^ tr [nia2 {Xt^XuXpXa + ^p^uip£.cr}] ■ 
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